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ABSTRACT 

Mission scenarios outside the Earth’s protective magnetic shield are being studied. Included are high usage assets 
in the near-Earth environment for casual trips, for research, and for commercial/operational platforms, in which 
career exposures will be multi-mission determined over the astronaut’s lifetime. The operational platforms will 
serve as launching points for deep space exploration missions, characterized by a single long-duration mission 
during the astronaut’s career. The exploration beyond these operational platforms will include missions to planets, 
asteroids, and planetary satellites. The interplanetary environment is evaluated using convective diffusion theory. 
Local environments for each celestial body are modeled by using results from the most recent targeted spacecraft, 
and integrated into the design environments. Design scenarios are then evaluated for these missions. The 
underlying assumptions in arriving at the model environments and their impact on mission exposures within various 
shield materials will be discussed. 


INTRODUCTION 

Human missions beyond LEO will first be confined to the Earth’s neighborhood and typically less than 100 
days duration. With growing infrastructure and operations, astronauts will become career space workers with space 
career durations approaching from one to two thousand days over a ten-year career (Tripathi et al. 2002). This near 
earth development will provide intrastructure and experience for ventures beyond Earth’s neighborhood first to 
Mars and beyond with missions of a 1000 days and more. The primary technology thrust to enable this 
development is new propulsion systems and space radiation protection (O’Keefe 2002). A triage of space-radiation 
health-maintenance strategies includes crew selection, radiation shielding, choice of propulsion technology, 
operational planning, and biological countermeasures (Cucinotta et al. 2001). The control of the radiation insult to 
specific tissues is the first line of defense against adverse health effects. The physical insult can be controlled in a 
limited way through material arrangement and material choices in spacecraft design (Wilson 2000). 

Fundamental to radiation-health risk assessment is the specification of the space environment and its 
transition through spacecraft materials to specific tissues. In the present paper we discuss the development of 
software to evaluate the radiation environment anywhere in the solar system and at any projected time including 
e val uation of material shielding properties. 

DEEP SPACE ENVIRONMENTS 

There are two basic sources of space radiation in the solar system: galactic cosmic rays (GCR) leaking into 
the heliosphere from local interstellar space and solar particle events (SPE) associated with acceleration within the 



transition region (shock) between a coronal mass ejection (CME) and the local interplanetary medium (Reames 
1999). The inward diffusion of GCR into the solar system is countered by convective currents due to the outward 
expansion of the solar wind (Parker 1965) resulting in an intensity reduction with decreasing helioradius. As to 
SPE, the CME volume expands with increasing helioradius reducing the shock intensity resulting in less intense 
SPE. Very large SPE intensities result within the shock region. When a large body is encountered, the local GCR 
or SPE environment is effectively shielded in the direction of the body, new particles (especially neutrons) are 
produced in interactions with the body atmosphere and surface, and if there is an associated magnetic field with a 
large dipole component then the neutron decay into charged particles (proton and electron) may be trapped leading 
to intense radiation belts with maximum intensity near the magnetic equator (Divine and Garrett 1983). 

GCR Environment 

A solution of the appropriate Fokker-Planck equation was found by Parker (1965) in which the inward diffusion of 
galactic cosmic rays is balanced by the outward convection by the solar wind. The density of cosmic ions within the 
solar system, assuming spherical symmetry, is then related to the external density as 

lKr,R) = Ho(R) exp [-/ r V(r’) dr 7 D(r’,R)] (1) 

where jl(r,R) is the ion density at radial distance r and magnetic rigidity R (particle momentum per unit charge and 
determines the trajectory radius of curvature in a magnetic field), Po(R) is the density in interstellar space, V(r) is 
the solar wind speed, and D(r,R) is the diffusion 
coefficient (Balasubramanyan et al. 1967). The 
wind velocity and diffusion coefficient depend on 
the solar activity usually measured by a 
subjectively defined number of sunspots (Wolfe 
sunspot number) observed in the solar surface (see 
discussion in Wilson et al. 1999a) and there is a 
phase shift between sunspot number and 
modulation as the wind generated at the solar 
surface diffuses into the modulation cavity which 
extends far out into the solar system. The relation 
of sunspot number to the cosmic ray induced 
neutron monitor count rate in Deep River, Canada 
is shown in Fig. 1 during some of the more recent 
solar activity cycles with projections to 2022. The 
development of the statistical model of sunspot 
number and its application is discussed elsewhere 
(Wilson et al. 1999a, 2002a), We note that the 
cycle is reasonably well determined soon after solar 
minimum but chaos enters late in the cycle and 
great uncertainty enters in cycle projections across 
the solar minimum (Wilson et al. 1999a). The 
uncertainty from this chaos is not included in the 
current projections. The inverse relation of solar 
activity and cosmic ray intensity is clearly seen in 
the figure. A simplified version of the Parker 
diffusion model was implemented by Badhwar et 
al. (1994) in which the solar wind is held constant 
at 400 km/s and the diffusion coefficient is taken as 
a function of lime and is correlated with the Mt. 

Washington neutron monitor count rate. The 
diffusion was found to be bimodal with unique 
dependence on the orientation of the solar magnetic 
dipole. Assuming an isotropic diffusion coefficient 
in which the diffusion coefficient generally 
increases with radial distance as D(r,R) = Do(R) r where s is on the order of 0 to 2, these assumptions lead to 

p(r,R) = jlo(R) exp {- V 0 (r 0 ** - r l s )/[(l - s)D 0 (R)]} (2) 



Fig. 1 Sunspot Number and Deep River Neutron Monitor 
countrate, actual and projected. 



Jupiter, and Saturn. 




where V 0 , r 0 , and D 0 (R) are the wind speed, size of the modulation cavity (50 to 100 AU), and diffusion coefficient 
at 1 AU respectively. Equation (2) is used to scale the modulated flux at 1 AU to arbitrary distance (Wilson et al. 
1999a). Modulation studies using various Pioneer, Voyager, and IMP spacecraft show variability of .v with solar 
cycle for some restricted energy ranges but the 
gross behavior for all energies above 70 MeV per 
amu is well represented by .v = 0.5 (Fuji and 
McDonald 1997). 

The solar modulation is less effective at 
solar minimum but the the radial dependence 
during solar maximum is apprediable. The GCR 
environments near Mars, Jupiter, and Saturn are 
shown in Fig. 2 during the 2002 maximum in the 
solar cycle. It is clear that the lower energy GCR 
are less depleted at the larger heliocentric radii 
during solar maximum. This radial dependence 
will place additional demands on protection near 
solar maximum at large helioradii. 

SPE Environment 

Solar cosmic rays or solar particle events 
(SPEs) were first observed as sudden short-term 
increases in the ground level ionization rate. The 
close correlation with solar flare events first 
identified them as originating in the solar surface 
plasma with eventual release into the solar system (Meyer et al. 1956). Thus it was assumed that the observation of 
solar surface phenomena would allow forecasting the possibility of such events (FAA 1975). Modern opinion 
considers the particle acceleration region not to be on the sun at all (Reames 1999). Rather large coronal mass 
ejections (CME) from active regions of the solar surface propagate into the interplanetary environment carrying 
along with them the local solar surface magnetic field frozen into the ejected mass (a good electrical conductor). 
There is a transition region (shock) between the normal sectored magnetic field structure of interplanetary space 
and the fields frozen into the ejected mass which forms a shock in which the interplanetary plasma is accelerated 
forming the solar particle event. The escape of the particles from the acceleration region is diffusion limited by 
particle-wave interaction with the interplanetary plasma so that a maximum intensity is implied. However, when 
the acceleration region passes the observation point the intensity is observed to increase by an order of magnitude 
to high levels in so-called shock events and no upper limit in intensity is known within the shock region (Reames 
1998). The SPE energies obtained in the acceleration process are related to the plasma density and velocity of 
propagation of the CME mass. To understand the SPE is then to understand the release of coronal mass and its 
propagation into interplanetary space relative to the observation point (Reames 1999). 

The SPE associated with Earth based ground level events (GLE) are intense at high energies and of greatest 
concern to human protection. The largest such event observed was the February 23, 1956 event with a GLE 
reaching 3600 percent over background. The second largest such event was the September 29, 1989 event reaching 
370 percent over background as shown in Fig. 3. It has been demonstrated that a design environment taken as 
4xSeptember 29, 1989 assures that the even larger February 23, 1956 event would not have produced a significant 
biological consequence (Wilson et al. 1999b). We therefore use the 4xSeptember 29, 1989 as a design 
environment. This is approximately the 99 percent environment of Xapsos et al (1999), The radial dependence is 
problematic. Reames has suggested that very little SPE radial dependence may occur from Earth to Mars. Others 
have assumed Gauss law dependence. We have assumed an inverse dependence on helioradius to 2 AU and Gauss 
law beyond. Clearly the radial dependence of the SPE for missons near Earth or Mars is of critical interest. 

INDUCED ENVIRONMENTS 

The GCR and SPE are given by the above-discussed formalism in space that is free from any near bodies. 
As one approaches a large object such as a planet or a satellite, the ambient radiation field is disturbed by 
interaction with the near body. This interaction has several possible effects. In all cases, the produced neutrons 
diffuse into space providing a neutron albedo (to a lesser extent other particles as well). If there is an atmosphere 
then the neutrons on the surface of the body arise in part from interactions with that atmosphere and in part from 
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Fig. 3 Deep River neutron monitor count rates during the solar 
particle events of October 19 and September 29 of 1989. 



those diffusing from below the local surface depending on the extent of the atmosphere. For example, it has been 
shown that neutrons on Mars surface with energies below 20 MeV are dominated by those produced in the bedrock 
or regolith (Wilson et al. 1999c). The water ice and C0 2 ice are less prolific in neutron production and becomes a 
feature of the neutron map as will be shown. If, in addition, the body has a dipole magnetic field then trapping 
regions can capture the neutron decay products allowing the buildup of intense electron and proton belts about the 
body as in the case of Earth and Jupiter. The Earth neighborhood has been discussed elsewhere and will not be 
further discussed herein (Wilson et al 2002a). 


Lunar Surface Environment 

The GCR and SPE in the near Earth environment 
impact the lunar surface producing secondary particles, 
which diffuse from the surface into the local environment. 
In addition the lunar mass shields the surface from 
directions intersecting the lunar terrain. On a smooth 
moon, the shielding effect is nearly a factor of two 
corresponding to all directions below the horizon. But, 
even local surface features can be used to protect human 
operations (Simonsen et al. 1990). In lunar orbit, the 
shadow cast by the lunar mass decreases and the neutron 
albedo decreases. The GCR environment on the lunar 
surface is shown in Fig. 4 at solar minimum and solar 
maximum. Similar results are shown in Fig. 5 for the solar 
particle event of September 29, 1989. Note, the particle 
intensity is completely dominated by the direct solar 
particles arriving from above the horizon. The difference 
between Figs. 4 and 5 is that the high-energy primary 
particles of the GCR are prolific in producing neutrons in 
comparison with the SPE event. 

Mars Surface Environment 

Galactic cosmic rays and solar energetic particles 
impact the martian atmosphere where interactions stop the 
low energy charged particles and modify the composition 
of the particle fields penetrating to the martian surface. The 
penetrating particle composition depends on the amount of 
overhead atmosphere. On impact with martian surface 
materials provide additional modifications in the field 
including the production of neutrons which diffuse from the 
surface materials back into the atmosphere (Wilson et al. 
1999c). The diffusive neutrons depend critically on the 
surface composition which in fact varies widely over the 
martian surface and is seasonally dependent. The 
computational geometry is shown in Fig. 6. 

The martian atmosphere is represented as C0 2 
distributed according to the Marsgram model and the 
regolith composition is given by the Viking data with 
elevation from the Mars mapping missions. The seasonal 
dependent water ice and C0 2 ice are also represented. The 
transport properties are represented by the HZETRN 
transport code and database with the low energy neutron 
transport handled by the computational procedures 
developed by Clowdsley et al. (2000). 

Over regolith, the surface neutrons are dominated 
by those diffusing from the surface materials. This 
component is reduced over the C0 2 ice cap and becomes 



Fig. 4 Lunar surface induced GCR environment 
during the 1977 solar miminum (full lines) and the 
1990 solar maximum (dashed lines). 



Fig. 5 The lunar surface environment during the 
September 1989 SPE. 



Fig. 6 Computational geometry for Mars surface 
environment (Simonsen et al. 1990). 









smaller still over water ice.,This dependence is shown clearly in Fig. 7 where calculations of the total surface 
neutron flux spectra is shown along with the contribution from the overlying atmosphere (forward component in the 
figure). This will have some advantage in that the surface materials can be surveyed if the diffuse neutrons can be 
detected in low Mars orbit. Such an orbital model is under development and the ability to infer surface composition 
will be tested. 

GCR attenuation in the 


martian atmosphere with a buildup of 
secondary neutrons in reaching the 
surface as noted in features in Fig. 8. 
Clearly the high elevations are marked 
by intense surface neutrons where 
GCR impacts are high due to the 
reduced overhead atmosphere. The 
impact with the surface materials 
continues attenuation with greatly 
enhanced neutron production in 
regolith or bedrock (Wilson et al. 
1999c), reduced in CO 2 ice, and 
especially in the water ice. Over 
regolith or bedrock, neutrons below 
20 MeV diffusing from the surface 
materials into the atmosphere fully 
dominate the surface neutron 
environment giving enhanced and 
depleted neutron regions depending 
on composition. Hence, there are 
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Fig. 7 Mars surface neutron environment with 16 g/cm2 C02 overhead 
and various surface material compositions. 


great differences in the surface 
neutron environment depending on the 
atmosphere overhead and the local 
surface materials. Polar surface 
materials change radically with season 
as seen in Fig. 8. 

The degree of spatial 
resolution using neutron detectors in 
Mars orbit depends on the 
transmission characteristics of the 
Mars atmosphere. We anticipate that 
sufficient spatial resolution will allow 
detection of these features in Mars 
orbits of the Odyssey spacecraft and 
such a model is being developed. 
Note that these methods are distinct 
from the usual approach of detecting 
surface water where the neutron 
thermalization process is the marker. 
In this our case, it is the Fermi motion 
and multiplicity of neutrons in nuclei 
of the surface materials that allow the 
>10 MeV neutrons to be a marker for 
surface properties. The atmospheric 
generated neutrons for the variable 



Fig. 8 Mars surface neutron environment above 10 MeV 
variation. Note that the color scale is not absolute and comp 
two hemispheres is only for qualitative features. 


atmosphere is a confusing factor in interpretation of such measurements but the current computational models 


should allow such an analysis. An ideal instrument is available on the Odyssey spacecraft in its 400 km high 


inclination orbit for such studies (Prof. I. Mitrofanov, personal communication). 




In addition to applications to Mars science, the models will allow the use of these models to assess the 
impact on human exposure and allow design considerations for all components in either surface operations, in lava 
tubes, or in Mars orbit. 


Jupiter Neighborhood 

The solar wind generated convective currents 
produce less effects on the GCR near Jupiter relative 
to Earth. In addition the jovian magnetic field traps 
particles with electron intensities to large distances 
from Jupiter (70 Rj). The electron flux spectrum near 
Callisto is shown in Fig. 9. Near a jovian moon, the 
intensities of all components are reduced except the 
induced neutron fields. We are developing an 
Anytime/ Anywhere software package for use in 
mission analysis and an example of a mission from the 
Earth/Lunar LI to Callisto and return is shown in Fig. 
10 for a fixed mass of four shield materials. The 
mission starts near an assumed solar maximum (2045 
AD) using the same projection model as in Fig. 1. 
The large exposure rates mid-mission are due to the 
jovian electron belts until arrival on the Callisto 
surface where there is shielding below the horizon. 
The continued increase of the environment on the 
return to LI is due to decreasing solar activity. This 
software will soon include the Mars models presented 
in the previous section. 

CONCLUSIONS 

The deep space environments developed for 
mission analysis use physical models to extrapolate 
the limited environmental data in both space and time. 
Induced fields are evaluated using high-speed 
transport models. Interesting dependence of local 
induced fields on atmosphere and ground composition 
are found which may be validated using Mars orbital 
data. The broader software being prepared will allow 
radiation exposure evaluation for arbitrary spectra in 
the future but limited to near Earth, near Jupiter, and 
interplanetary space for now. The near Mars 
environment will be added in the near future. 
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